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Abstract: Thermal and interfacial properties of a QGP droplet in 
a hadronic medium are computed using a statistical model of the 
system. The results indicate a weakly first order transition at a 
transition temperature ~ (160 ± 5) MeV. The interfacial surface 
tension is proportional to the cube of the transition temperature ir- 
respective of the the magnitude of the transition temperature. The 
velocity of sound in the QGP droplet is predicted to be in the range 
(0.27 ± 0.02) times the velocity of light in vacuum, and this value 
is seen to be independent of the value of the transition temperature 
as well as the model parameters. These predictions are in remark- 
able agreement with Lattice Simulation results and extant MIT Bag 
model predictions. 
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1 Introduction 



The formation of QGP droplet (fireball) is one of the most excit- 
ing possibility in the ultra relativistic heavy ion collision (URHIC) 
[1]. The physics of such an event is very complicated and to extract 
meaningful results from a rigorous use of QCD appropriate to this 
physical system is almost intractable though heroic efforts at lattice 
estimation of the problem has been going on for quite some time 
[2]. One way out is to replicate the approximation schemes which 
have served as theoretical tools in understanding equally compli- 
cated atomic and nuclear systems in atomic and nuclear physics in 
the context of QGP droplet formation. This approach lays no claim 
to rigour or ab-initio "understanding" of the phenomenon but lays 
the framework on which more rigorous structures may be built de- 
pending on the phenomenological success of the model as and when 
testable data emerges from ongoing experiments. 

Our intention in this paper is to reevaluate the free energy of 
a QGP droplet in a bulk hadronic (pionic) medium, again in the 
limit of vanishing chemical potential, but using a different semi- 
phenomenological model for the system. The MIT bag model is 
simplicity itself; it puts all quarks and gluons as free particles inside 
a bag and makes the impermeable bag as the agent of confinement 
by ascribing a set of boundary conditions for quarks and gluons. It 
is fine to use the MIT bag model to describe the hadrons as bags of 
quarks, antiquarks and gluons, but to extend the idea to represent 
the phase boundary between the QGP droplet and the bulk hadronic 
medium makes one a bit uneasy, though it is plausible. And, this is 
precisely the assumption made by the earlier authors [3] who have 
used the MIT bag model for the system of QGP droplet in a bulk 
hadronic medium. It is our intension to seek an alternative approach 
to the MIT bag model which does away with the assumption i.e; 
the confining bag of the hadrons has the same property as the in- 
terface separating the two phases, we propose an alternative model 
to represent the same physical situation. 

2 The template density of states of Thomas- 
Fermi and Bethe and their QGP variant 

In a very elegant and successful statistical model of atoms of large 
atomic numbers Thomas and Fermi [7] demonstrated the way to 
compute electronic density of states to very high order of accu- 
racy. The Thomas Fermi model of atom assumes the electrons to be 
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Fermi-Dirac gas confined within a localized region by the confining 
electrostatic potential V(r) of the central nucleus. The potential 
is assumed to be very slowly varying in the region with the aver- 
age thermal energy T (setting the Boltzmann constant to unity) is 
small compared to V(r) within the region and comparable to it near 
boundary. 

It is now straight forward [7] to compare the electronic density 
of states, assuming all states to be filled in a volume v 

Ne = PL* ^/SVT 2 (1) 

The maximum kinetic energy of the electron at any point in phase 
space should not exceed the electrostatic potential (confining) at 
that point and therefore 

p^ ai /2m = —V(k), when k is the phase point under consideration 
and V(k) is the momentum transform of the coordinate potential 
V(r). Therefore, the total density of states in phase space is given 
by 

p e {k)dk = [-2mV(k)f 2 u/3TT 2 (2) 



or, 



p e {k) = K2m) 3/2 /27r 2 ] [-V{k)\ 1 / 2 ■ 



dV{k) 



(3) 



dk 

In a modified statistical model, density of states for nucleons 
Bethe[7]gave a formula with a chosen potential as 

tt 1 / 2 r- 
= uk5/4al/4 exp{2Vak) (4) 

where 'a' is the density paprameter which satisfies 

for single particle level density. 

We can adapt the basic ideas of these models to the case of a QGP 
droplet, the electrons or nucleons get replaced by quarks which are 
also Fermions, and the minimum kinetic energy of the quarks at 
each point in phase space must exceed the confining/ de-confining 
potential at that point, since the QGP by definition is a deconfined 
gas of relativistic quarks and gluons as against the non-relativistic 
electron of the conventional Thomas-Fermi Model or the nucleon 
of the bethe model. Therefore, p min = [—V con f(k)] and p ma x = 
[ — Vconf (oo)] which represents a reference energy and can be set to 
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zero, remembering that we are dealing with a relativistic system 
and where l k' refers to the corresponding quark momenta in phase- 
space. So an expression similar to (3) holds for the quark density of 
states, with the replacement of V(k) with a suitable QCD induced 
phenomenological potential. The quark density of states therefore 
is 



p q (k) = (u/n 2 )[-V conf (k)}' 



dV conf (k) 



(6) 



dk 

In this adaptation of the Thomas-Fermi and the Bethe ideas, we 
only capture the spirit of the original idea for a system which is very 
different in detail. 

It is easy to observe that this is a typical first order phase transi- 
tion behaviour indicated by the Ramanathan et.al approach, thereby 
allowing direct comparision and the opportunity to use the calcula- 
tional advantages of each of the approaches to bolster the usability 
of the approximation schemes as a phenomenological tool in the 
analysis of the data as we shall see in the following. 

In the approximation schemes of the Ramanathan et. al [5], the 
relativistic density of states for the quarks and gluons is constructed 
adapting the procedures of the Thomas-Fermi construction of the 
electronic density of states for complex atoms and the Bethe density 
of states [7] for nucleons in the complex nuclei as templates. The 
expression for the density of states for the quarks and gluons (q, 
g)in this model is 



-W fc )) 2 (^^JJ gfl (7) 

where k is the relativistic four-momentum of the quarks and glu- 
ons, v is the volume of the fireball taken to be a constant in the first 
approximation andV con f is a suitable confining potential relevant to 
the current quarks and gluons in the QGP [5] given by 

= (V2*)7 M 9 2 (k)T 2 - ml/2k (8) 

where m is the mass of the quark which we take as zero for the 
up and down quarks and 150MeV for the strange quarks. The g(k) 
is the QCD running coupling constant given by 

g\k) = (4/3)(127r/27)(l/Ml + k 2 /A 2 )) (9) 

where A is the QCD scale taken to be 150 MeV. 
The model has a natural low energy cut-off at 
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k m in = ( 7M iV 1/3 T 2 A 2 /2) 1/4 (10) 

with 

N = (4/3)(12vr/27). 

The above cut off is equivalent to a deconfinement condition 
wherein the energy of the quark and gluon at each point in phase 
space has to exceed the value of the confining potential at that point. 

The freee energy of the respective case i (quarks, gluons, interface 
etc.) for Fermions and Bosons (upper sign or lower sign) can be 
computed using the following expression 

Fi = TT 9l J dk Pi {k) ln(l ± e -(V^)/ r ) (11) 
with the surface free-energy given by a modified Weyl [8] expression: 

Fsurface = -^R^T 3 ^ (12) 

where the hydrodynamical flow parameter for the surface is : 

7 = y/2 x y/(l/7 g )* + (l/ 7g ) 2 , (13) 

For the pion which for simplicity represents the hadronic medium 
in which the fireball resides, the free energy is : 

F w = (3 T/2tt 2 ) (-v) r k 2 dk ln(l - e -y/<+ k2 / T ) (14) 
Jo 

With these ingredients we can compute the free-energy change 
with respect to both the droplet radius and temperature to get a 
physical picture of the fireball formation, the nucleation rate gov- 
erning the droplet formation, the nature of the phase transition etc. 
This can be done over a whole range of flow-parameter values [5], 
We exhibit only the two most promising scenarios in fig.l and fig. 2. 
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Figure 1: F tota i at 7 5 = 87^, 7 9 = 1/6 for various temperatures. 
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Figure 2: F to tai at 7 S = 67^ , 7 9 = 1/6 for various temperatures. 
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3 Thermodynamical behaviour of the droplet and 
the nature of phase transition 

Standard thermodynamics gives the following relations: 

Entropy S = -(dF/dT) (15) 



Specific heat at constant volume Cy = T(dS / dT) v (16) 

Sound — Velocity C$ = S/Cy (17) 

After plugging the total free energy in the previous section into 
these expressions, we evaluate these quantities for the two promising 
scenarios of our model, the result of which are displayed in the 
following figs. 3 to 8. 
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Figure 3: Variation of S with temperature T at 7 9 = 67^ , j q — 1/6. 
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Figure 4: Variation of S with temperature T at 7 9 = 87^ , j q — 1/6. 
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Figure 5: Variation of specific heat Cy with temperature T at 7 5 = 

67 9 , lq = 1/6. 
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Figure 7: Variation of velocity of sound squared C| with temperature 
T at 7 9 = 6-fg , 7 9 = 1/6. 



13 




T (MeV) 



Figure 8: Variation of velocity of sound squared C| with temperature 
T at 7 9 = 87, , 7 9 = 1/6. 
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4 Surface tension of the interface 

The nucleation process is driven by statistical fluctuations being 
determined by the critical free energy difference between two phases. 
The Kapusta et. al model [4] uses the liquid drop model expansion 
for this , as given by 



AF 



47T 



R 3 [P had {T,fi B )-P qtg (T,fi B )} 
3 



3 

+4:% R 2 a + T crit T In 



1 + (^r)R 3 s q , g 



The first term represents the volume contribution, the second 
term is the surface contribution where a is the surface tension, and 
the last term is the so called shape contribution . The shape con- 
tribution is an entropy term on account of fluctuations in droplet 
shape which we may ignore in the lowest order approximation. The 
critical radius R c can be obtained by minimising (18) with respect 
to the droplet radius R , which in the Linde approximation [4] is, 



2(7 3AF C 



(19) 



T 


AF C 


R c 


a 


(7 

Ti 


(MeV) 


(MeV) 


(fm) 


(MeV/fm 2 ) 


150 


332.203 


3.475 


6.568 


0.078 


160 


382.359 


3.385 


7.966 


0.078 


170 


433.037 


3.285 


9.580 


0.078 


190 


532.219 


3.085 


13.35 


0.078 


210 


623.349 


2.875 


18.004 


0.078 


230 


702.041 


2.655 


23.776 


0.078 


250 


766.041 


2.455 


30.343 


0.078 



Table-1 for Surface Tension of QGP droplet at 7 g = 87^,7^ = 1/6. 
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160 


1197 


5.965 


8.031 


0.078 


170 


1494 


6.085 


9.633 


0.078 


190 


2216 


6.275 


13.435 


0.078 


210 


3088 


6.375 


18.140 


0.078 


230 


4059 


6.375 


23.844 


0.078 


250 


5052 


6.275 


30.630 


0.078 



Table-2 for Surface Tension of QGP droplet at 7 5 = 67^,7^ = 1/6. 

The tables 1 and 2 list the surface tension computed in the Ra- 
manathan et.al model which can be used in the dynamical models[4], 
thus reducing the arbitrariness in the latter models to this extent, 
thus enabling us to use the parameter extracted from a static sit- 
uation to make perdictions about the dynamical growth process of 
fireball formation. 

5 QGP droplet formation rate 

In this section we discuss the strength and weakness of two approx- 
imation schemes which seem promising in their usefulness in the 
above context and also explore their parametric interrelations which 
may be useful in phenomenological applications. The two approxi- 
mation schemes under consideration are the QCD oriented QGP- 
droplet formation model of Csernai and Kapusta and its later re- 
finements [4] and the simple statistical model of Ramanathan et. 
al [5] which again is essentially QCD oriented in that use is made 
of an effective QCD potential. In both these schemes it is possible 
to estimate the fireball radius, the transition temperatures for their 
formation, nucleation rate etc. in quantitative terms for compari- 
sion with raw data as and when it is available from ongoing URHIC 
experiments. 

As we shall see the two approaches can complement each other 
in the event of analysing the data on fireball formation especially 
with regard to nucleation of QGP droplets in a hadronic medium. 

Central to this comparision is the rate I to nucleate droplets of 
QGP in a hadronic gas per unit time per unit volume is given by [6] 

I = I e ~ AF ^ (20) 

Where Iq is the nucleation rate at vanishing change in free energy 
AF of the system due to the formation of a single critical size droplet 
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Figure 9: The free energy difference AF(R) between a hadronic phase 
with and without a quark-gluon plasma droplet [4] 

of plasma. In the Kapusta et. al [3] formulation it is not possible to 
estimate the whole prefactor in terms of other possible measurable 
parameters of the droplet. There is no way to estimate the value 
of the crucial input a within the Kapusta et. al approach and the 
related approaches, though a value of 50 MeV/ fm 2 is assumed. The 
critical free energy in this model varies as in fig.9 [4]. 

As could be seen the general nature of the curves given by figs.l, 
2 and 9, though using two totally different approches leads to the 
conclusion that both predict a first order phase transition for the 
fireball production process. As a derivative computation we can 
also compute the nucleation rate leading to the droplet formation, 
which for the two scenarios is illustrated by figs. 10 and 11. 
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Figure 10: The fireball formation rate with temperatures at 7 9 = 67^, 
l q = 1/6 • 
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As could easily be observed from figs 10 and 11, for both the sets 
of parameter values, the droplet formation rate grows exponentially 
in the vicinity of transition temperature range of 150 to 170 Mev 
as expected from lattice simulations and other model calculations. 
The critical free-energy and radius at the maximum of the curves 
in figs. 1 and 2 allow us to compute the surface tension used in 
the Kapusta et.al and the derivative models [4] as already shown in 
section 4. 

6 Conclusion 

The graphs (Figs. 3 to 6) clearly indicate that the model pre- 
dicts a weakly first order transition at a temperature in the range 
(160 ± 5) MeV, which seems to be consisitent with current ex- 
pectations of QGP-Hadron phase transition [1]. This feature is also 
corroborated by the behaviour of the nucleation rate exhibited in 
Figs 10 and 11, where the exponential rise in the droplet formation 
rate in the vicinity of 160 MeV, is indicative of the occurence of 
phase transition at these tempeartures. 

The independence of the velocity of sound in the QGP system, Figs 

7 and 8, from both the values of the model flow parameters as well 
as the magnitude of the transition temperature is remarkable. The 
value of the velocity of sound is consistently predicted to be of the 
order of (0.27 ± 0.02) times the velocity of light in vacuum. The 
tables 1 and 2 list the surface tension computed in the Ramanathan 
et.al model which can be used in the dynamical models[4], thus re- 
ducing the arbitrariness in the latter models to this extent, thus 
enabling us to use the parameter extracted from a static situation 
to make perdictions about the dynamical growth process of fireball 
formation. 

In both tables 1 and 2 the surface tension is seen to increase with 
the temperature of the fireball, which is a beautiful demonstration 
of a QCD effect. As the temperature of the QGP droplet increases 
the shear forces on the fireball surface will also increase tending to 
tear the surface quarks apart, consequently, bringing into play the 
confining property of the QCD forces manifesting itself in increased 
surface tension, which is exactly what the calculations show.Another 
striking feature of the result is the independence of the QGP droplet 
surface tension a variations in the flow parameters of the model and 
it varies with only the temperature, in the lowest order approxima- 
tion we have employed. 

The constancy of the ratio indicates a cubic crtitical temper- 
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ature dependence of the surface tension of the interfacial separation 
between the two phases. This is in striking conformity with the 
results of Lattice QCD simulation [9]. 

Finally it is readily seen from Tables 1 and 2, that at a transition 
temperature of 156 MeV, the surface tension is predicted to be of 
the order of 67 MeV/ fm 2 , a value in toatl agreement with the MIT 
bag model prediction of Abbas et.al [10]. 
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